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A single-electron transistor (SET) in a magnetic field irradiated with microwaves is studied the- 
oretically in non-equilibrium Kondo regime. The two fold effect of frequency-f2-microwaves is con- 
sidered as follows: the oscillations in the voltage with frequency fl and in the coupling parameters 
with frequency fl/p {p £ N). We describe the system by the Kondo model at a specific point in 
the Toulouse limit. A non-perturbative technique is proposed, namely, the non-equilibrium Green's 
functions and physical observables are averaged over a period of 2np/Q. When the microwave irradi- 
ation is considered affecting only the voltage, one sees the Kondo satellites as stated in the previous 
studies. Moreover, the features of the differential conductance and magnetic susceptibility of a SET 
become richer when the Kondo couplings are considered oscillating on time. We obtain the satellite 
peak splitting. It explains the possibilities one can find in experimental results that the distance 
between peaks, which appear in the differential conductance - magnetic amplitude characteristics 
G (H) or in the differential conductance - dc voltage characteristics G(Vdc), can be smaller than 

hn. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.40.Gk 



I. INTRODUCTION 

The Kondo effect, stemming from a macroscopic quan- 
tum coherent coupling between a localized magnetic mo- 
ment and a Fermi sea of electrons, has been used to ex- 
plain successfully many extraordinary properties of bulk 
dilute magnetic alloy systems P,[H for a long time. Based 
on the developments in mesoscopic physics and nano- 
technology, it has been proposed to study the systems 
of reduced dimensions and the single impurity problems. 
The latter is also represented by the very rapidly devel- 
oping fields of quantum dots (QDs), where Kondo-like 
impurities are manufactured [3|-|5|. The local magnetic 
moment is screened by hybridization with the delocal- 
ized electron-spins, leading to the formation of a bound 
spin-singlet state. A many-body resonance in the den- 
sity of states (DOS) appears at the Fermi energy, which 
strongly influences the conductance of the system at low 
temperatures. Above the Kondo temperature, thermal 
fluctuations destroy the coherence, resulting in a non- 
monotonic dependence of the conductance on tempera- 
ture. Under the Kondo temperature, due to the Kondo 
correlations, the system becomes dynamic and can be 
controlled experimentally by tuning external parameters 
such as source-drain voltage, gate voltage, and magnetic 
field 0. 

Since the Kondo effect is a typical example of a strongly 
correlated system, it is of particular interest in the field 
of non-equilibrium phenomena [1, 01 ■ As our understand- 
ing of the Kondo effect for a time-independent bias has 
increased Ml the theoretical effort has turned to time- 
varying fields [sHlJI ■ The ac field produces satellite peaks 
in the dependence of the differential conductance on the 
dc bias, as a consequence of the fact that the microwave 
field gives rise to effective changes in the DOS versus 
energy which are dramatically illustrated in the tunnel- 



ing current [151 . However, this interesting result could 
not be observed experimentally until 2004 by Kastner's 
group (l6| . They measured the differential conductance 
of a SET irradiated with high frequency microwaves. The 
results show that when the photon energy is somewhat 
greater than the Kondo temperature and the microwave 
voltage is carefully chosen, single-photon satellite peaks 
in the differential conductance can be observed. Very 
recently, Kogan's group has investigated the differential 
conductance as a function of magnetic field for different 
high frequencies of microwaves |17| . The non-trivial dy- 
namics that emerges when external parameters compete 
with the Kondo correlations has attracted our attention 
to study the Kondo effect in a SET irradiated with mi- 
crowaves. Due to the complexity of out-of-equilibrium 
many-body systems, it is instructive to investigate those 
cases where the exact non-perturbative solutions are ac- 
cessible. One particular interesting case is the Kondo 
model in the Toulouse limit, which has been solved ex- 
actly in equilibrium [l8| as well as out of equilibrium 
(SH theory) [l^. A question arising from the SH theory 
|10| . which is related to the experiment, is what happens 
if the Kondo coupling parameters are also affected by 
microwave irradiation. We can have different scenarios 
in different experiments: i, only the voltage oscillates; 
ii, only the Kondo couplings oscillate; iii, both of them 
oscillate [l^. Moreover, the multi-effect of microwave ir- 
radiation can induce the fact that the voltage oscillates 
with frequency fl, the Kondo couplings oscillate with fre- 
quency fii = ri/p, p e N in an appropriate experimental 
set-up. We predict that p is a frequency-dependent sys- 
tematic parameter. 

In this paper we investigate the response of a SET 
to a simultaneous application of a magnetic field and 
a high frequency microwave irradiation. We consider 
Kondo couplings oscillating with a frequency integer 
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times smaller than the frequency of microwaves. We 
solve the problem non-perturbatively by averaging non- 
equilibrium Green's functions in the period of the slow os- 
cillation. The differential conductance and the magnetic 
susceptibility are computed. In these results, the Kondo 
satellites appear due to the oscillation of the source-drain 
voltage. These satellites are split when the oscillation of 
Kondo couplings is taken into account. The peak split- 
ting distance depends on the difference between the in- 
put frequency and the frequency of the Kondo couplings, 
which is assumed to be a natural number p. 

This paper is organized as follows. In section HH the 
time-dependent Kondo model in the Toulouse limit with 
time dependent Kondo couplings J"^(t) is introduced 
and mapped onto a quadratic effective Hamiltonian. The 
averaged non-equilibrium Green's function method used 
to calculate the differential conductance and the mag- 
netic susceptibility is presented in section IIIII The re- 
sults for the differential conductance and the magnetic 
susceptibility are presented in section IIVI We conclude 
our results in section FVl 



II. TIME DEPENDENT KONDO MODEL 

A. Model and proposal for time-dependent Kondo 
couplings 

The system we study is a QD attached to two leads 
by high resistance junctions, so that the charge of the 
dot is quantized. The Kondo effect emerges in a QD 
occupied by an odd number of electrons at tempera- 
tures below the mean level spacing in the dot. Un- 
der such conditions, the system can be described by 
a single-impurity Anderson model with the level spac- 
ing playing the role of cutoff in it jl^l. Moreover, we 
consider the dot in the Kondo regime: both ionization 
and electron addition energy are much bigger than the 
tunneling rates: {Ed, U — E^} ^ ^l/b.^ s-i^d the ap- 
plied fields do not drive the dot out of this regime: 
{eVdc, eVdot, eVac, HsgzB} < {Ed, U - Ed}. Hence, the 
Anderson model can be mapped onto a two lead Kondo 
model through a time dependent Schrieffer- Wolff trans- 
formation 11, 13| . 

In our proposed Kondo model, the two non-interacting 
leads of spin-1/2 electrons, labeled right (R) and left (L), 
each subjects to a separate time-dependent voltage, in- 
teract via an exchange coupling with a spin-1/2 impu- 
rity moment ~^ placed in between them. The conduc- 
tion electrons in each lead, that couple to the impurity, 
are described by one dimensional fields -^J^^ (x), where 
a = L,R, and a \, are lead and spin indices, respec- 
tively. The interaction of conduction electrons in 
the leads and the impurity is local in space at a; = and 
involves the conduction-electron spin densities, which 
are written as = \ V-L (0) C^)a,a' ^P^' (0)- 



Thus, the time-dependent Kondo Hamiltonian describing 
a SET in a magnetic field and irradiated with microwaves 
is 
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V{t) 



E 



dx 



cx,^—L,RX—x,y,z 



(1) 



The system is driven out of equilibrium by applying 
a voltage bias across the junction, V (t) = Vdc + 
Vac cos {^t). The dc part Vdc fixes a chemical-potential 
difference between the two Fermi seas. The ac part 
Vac COS {fit) fluctuates these two chemical-potentials as 
an oscillation on time, allowing us to average the Green's 
functions, which will be shown in detail in the next 
section. We notice that we adopt atomic units with 
/i = fcs = e = /^B(7i = 1 in this paper. 

In general, due to the microwave irradiation, the cou- 
plings are split into two parts: the time- independent part 
and the time-dependent part, which oscillates on time, as 



Jx' it) = Jxo + Jxi COS {n,t + M 



(2) 



a, (3 — L, R imply the left and right leads. The oscilla- 
tion frequency 51i of couplings is different from the mi- 
crowave frequency f2. In order to enter the exact solvable 
situation, we assume [l^, iM 21 1 



(3) 



If the SET is in equilibrium or if a dc voltage is applied, 
the coupling parameters are constant and equal to J^q. 
In equilibrium case, from the scaling equations to lowest 
order in the couplings, we find that the constraints in 
equation ([3]) remain stable upon scaling if either 
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rLL _ jRR 
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(4) 



(5) 



As stated in Ref. [21|, these two cases correspond to 
two distinct isotropic two-channel Kondo model. We no- 
tice that in the constraint the two leads are decou- 
pled, the channels are just the right and left leads, which 
carry no current. When the system is irradiated with 
microwaves, both the voltage bias and the perpendicular 
couplings are proposed oscillating on time. The model in 
this paper is anisotropic in real space. It cannot be ob- 
tained by the direct mapping from the Anderson model 
to the Kondo model. However, the experimentalists can 
apply the relations between the Kondo couplings and the 
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tunneling rates, which are stated in Refs. [ll, 13 1, in or- 
der to compare their experimental results to the results 
in our theory. 

For simplicity, we propose the time dependent perpen- 
dicular coupling parameters as 



= Jt'\t) = jf^cos(r!it) , 



jf^(t) = jf (5' + jf ^ cos (rjit) 



Jf^(t) = J^^ + Ji"^ cos (OiO . 



tLL 



RR , 



(6) 



This proposal can be justified qualitatively as follows. 
The oscillation of Jit) comes from the oscillations of volt- 
ages such as source-drain voltage and gate voltage. The 
voltage biases oscillate with the same frequency fi, as the 
frequency of microwaves. The Kondo couplings thus os- 
cillate with frequency Vt. However, experimentalists can 
think about some experimental set-ups, in which one can 
control the frequencies of different voltages or the phases 
between voltage biases oscillate on time. The goal of this 
work is to have the Kondo couplings Jit) oscillating on 
time with frequency fii ^ VL. For simplicity, we assume 
either ili = pQ, or Oi = il/p, p G N. Since the VLi = VL/p 
assumption is predicted to give the results of satellite 
splitting, it is chosen to be studied first. The Vli = pJ7 
assumption will be represented later. 

Besides, one can investigate the time dependence of 
Kondo couplings in a general form as shown in equation 
However, the goal of this paper is to investigate the 
effects of the microwave irradiation. We thus consider the 
Kondo couplings as shown in equations Since we con- 
sider the time average of observables, we can neglect the 
phase difference among the coupling parameters. In gen- 
eral, microwave irradiation changes the transverse cou- 
plings Jj'^(t), thus changes the channel isotropy. We 
find that a QD either in a magnetic field [2^, or irradi- 
ated with microwaves can be driven out of the isotropic 
two channel Kondo limit. 



B. Mapping onto a solvable model 

The procedure to arrive at a solvable model is given 
in the following. We first bosonize the Hamiltonian ([T]) 
with notice that we have considered Klein factors and 
assumed the length was infinity: 



V'a<T(a;) 



(7) 



where ^aa{x) = a/tF J_dx'Iiaa{x') - (l)aa{x) 

4>aa{x) are Bose fields and Tiaa{x) are their conjugate 
momenta, satisfying commutation relations (23| . The 
next step is introducing new bosonic fields: charge 
4>c{x), pseudo-spin 4)s{x), flavor (f>f{x), and pseudo- 
flavor (psfix) as (jjcix) = [J2a,a'l^aa) /2, (t)s{x) = 



(Eo,.'^J>a.)/2, M^) = (Ea,. /2, 
4>sfi^x) = {j2a,a<a'^aa(l^aa) /2, and alsO for 

Hjy, <I>y, N,y, V = c, s, /, sf . We then perform 
the transformation of the Hamiltonian WHU^^with 
U = cxp [— ir^<i>s(0)]. Four more Klein factors Fy are 
introduced, which satisfy [Fy,Nyi] — 5,jyiFy and relate 
to the old ones as fI^Fl^ = F^F^f, F^^Fr^ = i^.Fj^, 
^Ri^Lt = FfF,, fI^Fli = FfFl The new impurity 
fermion is thus represented as d) = FgT^ ^ d = F^t^, 
= d^d — 1/2. We now re-fermionize these bosonic 
fields as 



1 



--F„ 



(8) 



to re-write the Hamiltionian in which all the charge, spin, 
fiavor, and spin flavor degrees of freedom are separated 

/oo 
dxtl'l{x)dxtpy{x) 

+ [H- (J, - 2tivf) : V'l(0)^'s(0) :] {Sd - 1/2) 

OO 

dx [Vdc + VacCOs{VLt)] 1pj{x)lpf{x) 



*t.(0)-*^(0)](d + d^) 

'^lfiO)-^,f{0)]{d + d^) 

*I/(0) + *./(0)k-d) , (9) 



where : Vl(0)V's(0) : means normal ordering with re- 
spect to the unperturbed -0,, Fermi sea, and we have 
defined Jt(i) = Jf^(<) = Jf ^ cos (r^it), J,(<) 
Wo +'^f(f)/2 + (^" + ./Dcos(f}iO/2, Ja{t) = 
[Jj^^ — J^^) cos (flit) /2. From our assumption in equa- 
tion 13) i Jz ~ 2t:vf, we are in the Emery-Kivelson 
line |l8| . the Hamiltonian ^ reduces to quadratic form, 
which is exactly solvable in equilibrium [18| and in the 




SH theory [10|, |21| . The charge and spin sectors are de- 



coupled from the local flavor and spin flavor, reduced to 
a collection of uncoupled harmonic oscillators. There- 
fore, they will be omitted from now on. Finally, by 
introducing Majorana fermions d ~ (d^ + d) / \f2 and 
h ~ [S — one arrives at the following Hamil- 

tonian 



Here 'H = ^ fk{t)c]f^Cfk + esfkclf^Csfk 



k 

-iHdb 



-Mfk - csk)a 



, iJsit) , t I \t , "^"(^) /- t ^- 



(10) 
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where efk{t) = 2ttv pk + e[Vdc + Vac cos (nt)], £sfk = 
2'KVFk. Here after we will work with Hamiltonian ([TO] 
to calculate physical observables. 



III. AVERAGE GREEN'S FUNCTIONS 

The Keldysh Green's function technique is well known 
as an efficient method to solve non-equilibrium problems 
(2^ . This technique has been applied to express the 
fully nonlinear, time dependent current through inter- 
acting and non-interacting resonant tunneling systems 
(ist - However, the complete result only obtained for the 
case level- width function (or tunneling amplitude) is time 
independent. Similarly, the Hamiltonian (llOp with con- 
stant couplings has been solved exactly [10|. Later, the 
steady state in the Kondo model in the Toulouse limit, in 
which the couplings are periodically switched on and off, 
has been investigated by analyzing exact analytical re- 
sults for the local spin dynamics at zero temperature [2^ . 
The exact solution for the time dependent Hamiltonian 
(jlOP still remains. In this section, we present our average 
Keldysh non-equilibrium Green's function method which 
gives good results in the high frequency regime. 



A. Non-interacting Green's functions of flavor 
fermions 

Because the first effect of microwave irradiation in- 
duces the oscillation in the voltage bias as shown in the 
Hamiltonian (jlOp. we first discuss the non-interacting 
Green's functions of the flavor fermion, which is defined 
as 

9fk{t,t') = -i(TK{cfk{t)c:'fi,{t')]) . (11) 

The flavor fermion is considered as in the two dimensional 
fermion sea. The retarded, advanced, and Keldysh com- 
ponents are 



g'^f {t, t') ^ Tie (±< T t') e-'2^^^K*-*') 

X g-i%i[sin(m)-sin(Ot')] 



9% 



{t, t') = i [2/ {2^VFk + Vdc) - 1] e-*2-.Kt-t') 

X g-«-%^[sin(m)-sin(Ot')] ^ j--[^2) 



where / is the Fermi distribution function in the flavor 
lead. While the chemical potential of the spin flavor 
lead is chosen as the reference, the chemical potential 
of the flavor lead oscillates on time with an amplitude 
Vac around the fixed chemical potential Vdc- Using the 
identity exp [x (a — 1/a) /2] = X]^-oo where 
J„ (x) are the integer Bessel functions of the first kind. 



and changing variables 

T = t-t' , 
t + t' 



T 



(13) 



we find 



m,n— — 00 



X e 



-i[ek+{n+7n)n/2]T -in{n-m)T 



(14) 



Because the fiavor fermion chemical potential oscillates 
on time with frequency 17, one can take the average of the 
non-interacting Green's function over time T ^ f2~^[23|. 
The Fourier transform regarding the time difference r of 
the average Green's function is 



9f [^) = Z^9fk (^) 



Q J ^-^ Lo — 2TrvFk — nQ 

m,n— — 00 ^ k 



The DOS of the flavor fermion is lli 



(15) 



Pf 



2 I Vac 



p{lo — nil) 



1 



1 



Vac 

Q, I 2'hvf 2'hvf 



(16) 



We have the second line because we have applied the 
linearization for the energy spectrum around the Fermi 
level [23. The retarded and advanced components thus 
remain as those in equilibrium 



9/ (w) = T 



2VF 



(17) 



However, the effect of ac voltage Vac {t) modifies the 
Keldysh component gj^ {t,t') as 



VF 



J2j?.{^][2f{u^-nn + Vdc)-l] , (18) 



with 9^ (uj) is the average of the Keldysh component in 
the Fourier space corresponding the time difference r. 

We then calculate the combined non-interacting 
Green's functions of the flavor fermion. We show here 
one example 



(19) 
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Its retarded, advanced, and Keldysh components are with 

mf^{uj) = 0, 



(a;) = (2^/.^)^J^(% 



X {/ - n^! + Vdc) - / (c^ + nn - Vd,)} ■ (20) 

About the spin flavor fermion, its chemical potential 
is chosen as the reference, and no oscillation voltage is 
applied on the spin flavor lead, it is considered in equi- 
librium. 

B. Interacting Majorana Green's functions: 

In order to compute observables, it is necessary to cal- 
culate advanced, retarded, and Keldysh components of 
Majorana Green's functions, which are defined as 

Go,p{t,t') = -i{TK{a{t)l3{t')}) , a,l3 = a,b. (21) 

We show here the results of the Green's functions 
G^^{t,t') and G^^^{t^t'), which are used to calculate the 
differential conductance and the magnetic susceptibility. 



1. Advanced Green's function Gaa{i,t'): 

We write the Green's functions in terms of the 
interaction-picture operators by invoking the S matrix, 
we have the advance Green's functions as follows 



dhJt{h)gCt{t-h)Gt{h,t'), 



+ 



AiravF J- 

/oo 
dhgt{t~h)Gt{ti,t') 
-OO 

/ dh [J? ih) + Jl [h)] gtit - h)Gt{ti,t') . 

J —OO 

(22) 



With the assumption of coupling parameters Ji (t) 
from PU)) . we have 

/OO 
dtMti-t)Gt{ti,t') 
-OO 

/OO 
dtiTs {h)&{ti~t)Gt{ti,t') , 
-OO 

/OO 
dhQ{ti~t)Gt{h,t') 
-OO 

/OO 
dtiTat{ti)e{ti-t)G^^{ti,t') , (23) 



r,, {t) = Tso + 2yr^cos (nit) + r, cos^ {n,t) , 

Tat (t) = Tat cos^ (nit) , (24) 
where we have defined Tt = ( J^'")'^ /Airavp, Ta = 



(Jf^ - Jf «) /mnavp, Tat ^ Tt + T,, r,,o = 

{Ji^ + J^^y /WiravF, r, = (Jf^ + jf«)Vl67rawj.. 
Equations (j23p induces the SH theory limit when all cou- 
plings are time independent. One can easily take Fourier 
transform of equations psp . then re-obtain the exact for- 
mulas of the advanced Majorana Green's functions, as 
shown in the formula (4.9) in Ref. plj . 

We need to solve the integral equations (|23|) and find 
the Green's function G^^{t, t'). We first change the equa- 
tions ([23|) to differential equations by taking the deriva- 
tive with respect to time t. We have 

[dt~T, {t)]GtM = HGtit,t') , 

[dt - r,, it)] Gt {t, t') + HGt (t, t') = -^5 (t t') . 

(25) 

We define the non-interacting Green's functions 
9%, (t,i')as[2l 

A 



[dt-T,t/, it)] g^i. it,t') = 6it-t') 



(26) 



so 



gtt{t,t')^-e{t' -t)cxp 



X exp 



^ [sin(2r2i<) - sin (20it') 



2 



{t t') 



X exp 



2vT^or^ 



[sin(r2it) - sin(17ii')] 



r 



[sin(2rJit) - sin(2r2it') 

4r2i 



X exp 

From equations (p5|) and ([26)) . we have 



(27) 



dti 



[g^t it, h)] " + H^g^ {t, h)\ Gt ih^t') = -^S [t t') . 

(28) 



2. Keldysh Green's function G^^{t,t'): 

In the same way as we did with the advanced Green's 
functions, we write the Keldysh Green's functions by in- 
voking the S-matrix as 

/OO 
dhHe{t~h)G't{h,t') 
-OO 

- / dtirs{ti)e{t-ti)Glf,{tut') + HTs{t,t') , (29) 
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/CO 
dUHe{t-t,)Gt[h,t') 
-oo 

oo 

dhTat{ti)Q{t~U)Gt{h^') + Tat{t,t') , (30) 



where we have defined the Tat {t,t'), Tg {t,t') as 

/•OO 

X [2/M - l]e— (*^-*^) / dts9t{t2,t3)Gtih,t') , (31) 



Jt ih) Jt (i2) 



/oo 

X [/ (w - nn + Vda) + f{LU + nn- Vdc) - 1] 



dtidt2- 



-e(t-ti) 



X / dc. [2/ (..) - 1] e— (*-*^)Gi(t2, i') . (32) 



The above equations induce the following differential 
equations 

[dt + r, (t)] Gtit, t') = HGtit. t') + HKs (t, t') , 

[a* + r,, {t)] G^Jt, t') + HG^,{t, t') = Kt {t, t') , (33) 



We realize that the non-interacting Green's functions 
g^tf^ (tyt') oscillate on time with the frequency fli = 
Q./p. We cannot apply any cut-off for the boundary of the 
integrals in the equations ([21]) and . It means we can- 
not solve these equations exactly. However, g^^l^ (t, t') 
can be considered oscillating around smooth averaging 

functions g^l^^ (r, r)with T :$> fl^^. It is discussed in 
detail in the following subsection. 



C. The Averages: 

In order to find the Majorana Green's functions 
G'^^{t,t') and G^^[t,t'), one needs to solve equations 
and ([55]) . In this subsection, we represent our aver- 
aging approximation method. 



1. Advanced Green's function Gaai<^) ■ 



The .g^j/^ {t^t') ii^ equations ([?7|) can be expressed as 
functions of the Bessel function products, for instance, 

..1(M') = -e(t'_t) &(-.^)^™(-^) 

(37) 



m.n— — oo 



X g%^(t-*')e'^^i("*-™*'). 



with Aat/s{t,t') = dtTat/s{t,t'). We define non- We change time variables such as r = < - f , T = 
n^^^^'o „K f-i (t + t') /2 . Because the Green's functions qC'^,(t,T) os- 



interacting Green's functions g^^^^ it^t') 

[dt + T,t/s{t)]g^t/At,t')-^S{t't'), (34) 



{t + t') /2 . Because the Green's functions T) os- 

cillate on time with period 27r51]^^, we average them over 
time T in a period of 2tt /^li. We then write it in Fourier 
space as 



5^ (t,<')-0(t-t')exp 

Tat 



{t t') 



X exp 



gfit,t') = e{t~t')exp 



sin(217ii) - sin(217it')] 



sO 



X exp 



2vI\oTs 



[sin(17it) — sin(r2it')] 



[sin {2nit) ~ sin (217it')] 



X exp 
We then have 

/oo 
dt:5f(t,ti)[Gfa(ti,t') + A.(ti,t')] , 



(35) 



dti 



— OO 
-1 



(36) 



i/at 



(38) 



Similarly, 



-2)i+2fe — * 



.2A/rsor.s 



f71 



m,n,A;— — oo 



cj + (m + 2fc) Oi - i (r^o + ^) 
We thus have 



g^.ico) +H^g^{uj) 



(39) 



(40) 
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2. Keldysh Green's function G^^iu)): 

In the same technique, we obtain the average non- 
interacting Green's functions 



2 I ^^at 



n— — oo 



r.- ' (41) 



+2fe t 



m.n.A;— — cxD 



Lo + (m + 2fc) Oi + i (Fso + ^) 



(42) 



We also average Kat/s {tjt') and obtain Aqj/^ (lu) as 



+f{uj +in- 1/p) n - Vdc)+ fioJ - (n - 1/p) n + Vdc) 
+f iuj+{n + 1/p) n - Vdc) - 2]} , (43) 



A, (u) = (co) Gt H {2r,o [2/ [u) ~ 1] 

+r, [/ {uj - n/p) + j{uo + n/p) - 1]} . (44) 

In this limit, G^(i, t') only depends on t = t — t'and its 
Fourier transform is 



Hg^iio) g^,{u;)AatH + As{io) 



1 + H^g^^ iuj)g^ (w) 

The meaning of our averaging method is that we have 
smoothed the fine oscillations of the Green's functions. 
It works well in the high frequency regime, namely, 
T ':$> . It offers the opportunity to compute the phys- 
ical observables non-perturbatively. Our results should 
be close to the exact solution while the results from the 
perturbation calculation cannot be. Any cut-off in the 
perturbation calculation will induce a loosing of some 
contributions in the physical observables. However, our 
averaging method cannot be applied for the adiabatic 
limit. The adiabatic regime should be considered sepa- 
rately, maybe as in Ref. [2§|. In the next section, we 
calculate the differential conductance and magnetic im- 
purity susceptibility by using this averaging method. 



IV. AVERAGE PHYSICAL OBSERVABLES 

In this section we will compute average differential con- 
ductance and magnetic susceptibility, and discuss the re- 
sults. We notice readers that the word "average" in this 
paper mostly refers our averaging method. 



A. Average charge current and average differential 
conductance 

We are first interested in calculating the time aver- 
aged charge current through the junction (Ic) because the 
time averaged differential conductance G = d{Ic) / dVdc is 
accessible experimentally We compute it by using 

the Keldysh non-equilibrium Green's function technique. 
The current at a time t can be written in the form 



Ic (t) = - 



eJtjt) 



(46) 



where we have defined the Green's function G^^.^^ (t, t) as 
Gfka{t,t') = ~t(TK{{c\,{t)+Cfk{t))a{t')}). 

The current is then expressed as Ic (t) = 
- ieJt{t)/8Tra)nJ^^ dtiJt{ti)mf {t,ti) G^ (h^t). 
As we apply the average condition to the Green's 
functions, the current at a time t means "average": 
IJt^^mf {t,h)G±{h,t)^mf {tM)Gt{ti,t) then 

Sir^avF J-oo J-oo ^ \^ 

X [/ (c^ + nfi - Vdc) -fiuj-nn + Vdc)] 
xe--(*-*^)Gi(ti,t) . (47) 

One should notice that we calculate the time averaged 
differential conductance, we therefore can apply our aver- 
age Green's function method, which has been discussed in 
section Inil The Majorana Green's functions G^^ {ti,t) is 
calculated averagely in the above section with its Fourier 
image is shown in equation (|40| . As discussed above, the 
current should oscillate on time with period 27ril^^. It is 
possible to take the average of the current over time t in 
a period of 2TTpfl~^. So the average current is 



— n 



2Trp/n 



Ic--7^ dtlcit) 



/oo ^ 
-oc- 



Vac 

n 



X {[/ iLj + in-l/p)n~Vdc) 
-f [lo - [n + 1/p) n + Vdc) 
+/ (io + (n + 1/p) n - Vdc) 
-f{co-in-l/p)n + Vdc)] Gtiuj)} 



(48) 



By taking the derivative of the time averaged current in 
equation (P5|) with respect to the dc voltage, we obtain 
the time averaged differential conductance. The conduc- 
tance and the magnetic susceptibility behaviors remain 
when the temperature is varied from absolute zero to a 
very small temperature. If the temperature is kept in- 
creasing much below the Kondo temperature, the con- 
ductance and the susceptibility behaviors are broadened 
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and then smeared. For convenient computation, we con- 
sider the time averaged differential conductance at the 
absolute zero 



77 



oo 

Sir ^ 



r 



Vac 



X [Gt[Vdc'{n-llp)n) 

+G^a{-Vdc + {n + l/p)n) 
Gt(yac~{n+l/p)n) 
+Gt{-Vdc + {n-l/p)n)] , 



(49) 



where the differential conductance is defined as G = 
2^ Jq'"^^^^ Ic (t) dt /dVdc- We wiU discuss the behav- 
ior of the average differential conductance as a function 
of magnetic field and source-drain voltage. 



B. Average impurity magnetization and 
susceptibility 

While the experimentalists are interested in measuring 
average differential conductance, the theorists are inter- 
ested in calculating magnetic susceptibility as well be- 
cause both observables provide direct information about 
the magnetic impurity of a SET. The local magnetic sus- 
ceptibility - the response to a magnetic field coupled to 
the impurity spin only - is the derivative of the magne- 
tization, which is calculated from the following formula 



M{H) = ^I d.Gt[u) 



(50) 



The Majorana Green's function (w) is calculated av- 
eragely in the above section with its Fourier image shown 
in equation (|45p . The magnetization is understood as the 
"average" one and is expressed as 



M{H) 



47r 



duj- 



{^iB9^) H f°° 9s (w) [9at (^)Aat (t^)+A, (w) 



(51) 

From formula ((5T|) . when the system is applied either a 
constant voltage bias V^d = Vdc or a time dependent one 
Vsd = Vdc + Vac COS {Qt) and the Kondo couplings are 
time independent, again, we re-obtain the results in the 
SH theory ^ 



21| 



We compute the impurity susceptibility as x (H) ~ 
fJ-BgidnM (H) and discuss its behavior in the following 
subsection. 



C. Results and discussions 

What we show below is just a sample of some plots to 
jive a flavor of our results. In all plots, the horizontal 
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magnetic field fiBdi^ 

Figure 1: (Color online) G (H) are plotted for Vac ~ 4, Fat = 
0.3, Ts = 0.05, Ft = 0.25, Vdc = 0, and ftfi = 2 in the up- 
per panel and h^l — 1.5 in the lower one. The dotted black 
curve is without microwaves; Kondo peak at HBQiB = 0. 
The dashed black curve is with ac voltage, constant cou- 
plings: J"^ = + J"^; Kondo satellites at fiBQiB — nMl; 
the highest peak's position depends on eVac/hO,. When 
both voltage and couplings are oscillating, peaks occur at 
fiBQiB = (n ± l/p)hQ; the solid red, blue, green curves are 
corresponding for the p = 1, 2, 3. The peak splitting due to 
time-dependent couplings can be seen clearly at the highest 
peak. 



axis is the magnetic energy H = fiB^iB, B is the ampli- 
tude of applied magnetic field, the vertical axis is con- 
ductance G / {e^ /h) or magnetic susceptibility x/idif^B)'^- 
We choose F^o = 1 as the energy scale and it can be 
considered Kondo temperature. It means the other pa- 
rameter values are chosen in comparing with Fso- The 
results obtained at very small temperature and at abso- 
lute zero are similar. For efficient numerical calculation 
we choose to work with conductance at absolute zero and 
susceptibility at very small temperature (T = 0.05). 

In figures [1] and |3l conductance G and susceptibility 
X are plotted as functions of magnetic field H for differ- 
ent situations: no microwave (dotted black curves), time 
dependent voltage (dashed black curves), and both time 
dependent voltage and couplings with different values of 
p (red curves for p = 1, blue curves for p ~ 2, and green 
curves for p = 3). When we compare the conductance 
with and without the effect of microwaves, we find that 
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microwave irradiation reduces strongly the conductance. 
Because we consider a SET in a special situation of cou- 
plings, we cannot compare our results with the case where 
we only consider oscillation in voltage Vac {t), no time os- 
cillation applied on couplings as we discussed in section 
|TT](no current in this situation in equilibrium). However, 
we can compare our result with the case in which we 
only consider oscillation in voltage Vac (t) and constant 



— r„, = 0.3, r,, = 0, r, = 0.25 

— r„, = 0.3, r, = 0.2, r, = 0.25 



couplings such as J"'^^ = J' 



±0 



7*^/3 7 J 

Jj_^ so Jt, JsO 



The side-band effect of microwaves with the peaks ap- 
pearing at ^B9iB = eVdc i nMl (n £ Z) due to the 
oscillation of source- drain voltage Vac if) [1, [IMll- As 
it is mentioned in Ref. [3| the photon-induced satel- 
lites appear quite delicate, we deeply investigate the ef- 
fect of ratio eVac/^^- The side-band peaks only ap- 
pear when eVac/fi^ is big enough compare to 1 (but 



eVa, 



hV,, as well as ftfJ ^ ksTx and hil > ksTx 



in order to avoid decoherence) , and this ratio determines 
the height of side-band peaks. For instance, in the up- 
per panel of figure [T] plotted for 1 < eVac/fi^ < 2, we 
see the highest peak at HBdiB = h£l besides the trivial 
main peak at nsgiB = 0. The strength of microwaves 
is small compared to the strength of the magnetic field, 
which reduces strongly Kondo conductance. When we 
increase the oscillation amplitude, microwave effect be- 
comes stronger than the effect of the magnetic field, it 
reduces strongly the peak at fJ.B9iB = and increases 
the satellite peak at ^BdiB = 2hil. as we see in the 
lower panel for eVac/h^ > 2. One can say the high- 
est satellite peak position depends on the ratio eVac/hV,. 
However, it is impossible to see the highest satellite peak 
at fJ-sgiB = 3hfl due to the strong suppression by both 
magnetic field and voltage. This result suggests to exper- 
imentalists an efficient method to calibrate the oscillation 
voltage amplitude Vac- 

We now discuss the second effect of microwaves. In 
fact, the oscillation of Kondo couplings comes from the 
oscillation of source-drain and gate voltages. It can be 
considered as an indirect effect. In figure [T] the dashed 
black curve is plotted for the constant couplings, the red, 
blue, and green curves are plotted for the time dependent 
couplings with p = 1, 2, 3 correspondingly. We find that 
when the coupling parameters oscillate with frequency 
Hi = fl/p with p G N, each main peak splits into two 
peaks. Since the energy is conserved, when a peak is 
split into two peaks, the height of these two peaks relate 
to the height of original peak. The distance between an 
original peak and its split peak is dpeak = Ml /p. This is 
understood mathematically by investigating the Green's 
function Gaai^)- The difference between the coupling 
frequency and the applied frequency is described by an 
integer number p > 1. However, p cannot be too big be- 
cause the slow oscillation in couplings can be considered 
as constant compared to the fast oscillation of voltage. 
When p ~ 1, the couplings oscillate on time with the 
same frequency as microwaves, the satellite peaks split 




r„, = 0.25, r, = 0.05, r, = 0.25 
r„, = 0.35, r, = 0.05, r, = 0.25 




I 

345 
magnetic field /i^gjB 

Figure 2: (Color online) G (H) are plotted for Vac = 4, HQ = 
1.5, ,p = 3, Vdc = 0, and different set of energy scales Fi. 
In the top panel. Tat = 0.3, Ft — 0.25 is fixed, Fs is varied: 
Fs = for the black curve, Fs — 0.2 for the red curve. In 
the middle panel. Fa = 0.05, Ft = 0.25 is fixed, Fat is varied; 
Fat = 0.25 for the black curve, Fs = 0.35 for the red curve. 
In the bottom panel. Fat ~ 0.3, Fs = is fixed. Ft is varied; 
Ft = 0.2 for the black curve. Ft — 0.3 for the red curve. 
Notice conductance amplitude strongly depends on Ft . 



into two peaks with the distance dpeak = fi^ equal to 
the satellite peak distance. So for p = 1, the peaks re- 
main but their height changes drastically due to the re- 
contribution of energy through the peak splitting. This 
peak splitting effect can be seen clearly if p = 2, or p = 3. 
As we have explained the number p can also be bigger 
than 1 due to the delay of phase between the source-drain 
voltage and gate voltage. The peak splitting can be seen 
clearly at the highest satellite peak. When we increase 
p the splitting distance dpeak decreases, the splitting be- 
comes weak then approaches the limit in which the cou- 
plings are constant. This result explains well a possibility 
in experimental result in which one sees peaks appearing 
at a distance iisgiB < Ml in the regime of weak magnetic 
field. We predict p is a frequency dependent systematic 
parameter. It should be checked carefully in experimen- 
tation. For the case p = 1 , we can check the height of the 
satellite peaks to confirm the two fold effect of microwave 
on voltage and couplings. 

Besides, from figure [21 we find the main characteris- 
tic of G {H) lightly depends on scales F^, (i = at,s,t). 
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Figure 3: (Color online) G (H) and x (H) are plotted in the 
upper and lower panels correspondingly for hQ — 1.35, eVac = 
2, Vdc = 5, Tat = 0.3, Ts = 0.1, Ft = 0.25. The dotted black 
curve is without microwaves; Kendo peak at fiBQiB = eVdc- 
The dashed black curve is with ac voltage, constant couplings: 
jap ^ _|_ ja^t. p-QjjjjQ satellites at ^ggiB = eVdc ± nftfi. 
When both voltage and couplings are oscillating, peaks occur 
at HBQiB — eVdc ± (n ± 1/p) hQ,\ the solid red, blue, green 
curves are corresponding for the p = 1, 2, 3. 

Ti are chosen much smaller than Fgo because the mi- 
crowave effects on the couplings should be small com- 
pared to the existing couplings. We vary each en- 
ergy scale Fs, Fat, Fj in the top, middle, and bottom 
panel correspondingly while other parameters are fixed. 
Ft ~ {J±_^^ concerns the coupling between two leads and 
it controls the spin-flip tunneling of electrons through 
QD. Fs ^ \^J±^ + ./j^^] concerns the sum of separated 
couplings in each lead and is small compared to Fso. 
(jj;^ — Jj^^) contributes to Vat- Both energy scales 
Vat and Vg affect slightly the G (H) behavior. The red 
curve in the bottom panel is plotted for the case in which 
modulation of microwave only affects the transverse cou- 
plings, so Fat = 0.3, Vs = 0, Ft = 0.3. The amplitude 
of conductance strongly decreases when we decrease Ft. 
It vanishes when Ft = 0. Thus, any small effect of mi- 
crowave irradiation drives the system away from the un- 
coupled isotropic two-channel Kondo situation. 

In figure [3l the magnetic susceptibility is plotted as a 
function of magnetic field at temperature T = 0.05. As 
previously stated from G (H) characteristic, dc voltage 



C5 



0.08 
0.06 


no microwave 
constant coup 

— p=l 
p = 2 

— p = 3 


ings 

( ' 
1 • 
1 ' 
1 ' 








1 


0.04 


i\ 
ih\ 


I 1 
I ' 
1 ' 
( ' ; 
.1 ' ; 






. ; 

■ ' lA 


•\ — 

M 

1*1 

'11* 


0.02 


'iV 
'IP 
. ( s 
— A ./ 11/ 


ml fV 






\t\ 11 


' 1 
Ml 

M 1 r A - 














-4 -2 2 4 
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Figure 4: (Color online) G (Vdc)is plotted for HQ = 1.3, eVac = 
3, Fat = 0.3, = 0.05, Ft = 0.2, H ^ 0. The dotted 
black curve is without microwaves; Kondo peak at eVdc ~ 0. 
The dashed black curve is with ac voltage, constant cou- 
plings; Kondo satellites at eVdc ~ inh^l. When both volt- 
age and couplings are oscillating, peaks occur at eVdc = 
± (n ± 1/p) fiJl; the solid red, blue, green curves are corre- 
sponding for the p = 1, 2, 3. 

Vdc also splits the satellite peaks Q when it is smaller 
than the ac voltage Vac- It may cause confusion to recog- 
nize the peak splitting due to the oscillations of couplings 
in experiment if experimentalists cannot avoid a small 
contribution of Vdc- The fact that Vdc splits peaks by a 
distance ±eVdc is understood in non-equilibrium Kondo 
physics of a QD without microwave irradiation (2l| . Sig- 
nificantly, if Vdc is bigger than Vac, we see the peak split- 
ting due to the oscillations of couplings. If Vdc = 0, we 
see the satellite peaks dominate at HBQiB = 0. The other 
peaks at nsgiB = kMl (fc = ±1, ±2, . . .) are very weak. 
To study the satellite peak splitting due to coupling os- 
cillation, we choose Vdc > Vac- We find all susceptibility- 
magnetic field and conductance-magnetic field character- 
istics are the same. It is evident that the relation between 
susceptibility and conductance shown in formula (8.9) of 
Ref. [21| can be generalized when the system is irradiated 
by microwaves. 

The differential conductance G and magnetic suscep- 
tibility X are also plotted as functions of dc voltage Vdc 
in figures 2] and [S] In order to investigate microwave 
irradiation effect, we consider the case in which no mag- 
netic field is applied. Again, we find the two-fold peak 
splitting (Kondo satellite splitting) due to both voltage 
and Kondo coupling oscillations. There is more chal- 
lenge for experimentalists if we are in the cases p = 1 
and p = 2. The distance between peaks is still h^l. For 
p ~ I, the positions of peaks are the same as those in the 
case when only voltage oscillates, but the peak heights 
are re-contributed. For p = 2, both the positions and 
heights of peaks are re-contributed, satellite peaks are 
shifted by Ml/2. For p > 3, our interesting results show 
that one has a chance to see two peaks, whose distance is 
much smaller than Ml. For instance, in the green curves 
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dc voltage eVfj^^, 

Figure 5: (Color online) x (Vdc) is plotted for Ml = 1.5, eVac ~ 
2.5, Tat = 0.3, r, = 0.1, Ft = 0.25, H = 0. The 
dotted black curve is without microwaves; Kondo peak at 
eVdc = 0. The dashed black curve is with ac voltage, con- 
stant couplings; Kondo satellites at eVdc ~ ±nftf2. When 
both voltage and couplings are oscillating, peaks occur at 
eVdc = i (n i 1/p) hQ; the solid red, blue, green curves are 
corresponding for the p = 1, 2, 3. 



in figures m and [SJ which are plotted for p ~ 3, one can 
find two peaks which are split from a satellite peak at a 
distance of 2hfl/3. Moreover, one can also find two peaks 
at a distance of fifi/S. They are the peaks, that are split 
from the consecutive satellite peaks. 



V. CONCLUSIONS 

In this paper, we have investigated the effect of mi- 
crowave irradiation on a SET in a magnetic field. The 
frequency of the microwaves is fl. The system is de- 
scribed by a non-equilibrium Kondo model in a specific 
point in the parameter space, where the Hamiltonian H 
can be diagonalized. Besides the oscillation in the source- 
drain voltage, the oscillation in the Kondo couplings is in- 
vestigated for a specific situation. The Kondo couplings 
J±^{t) relate to the source-drain voltage and the gate 
voltage, which oscillate with frequency 51. The couplings 
thus oscillate on time with the same frequency or its har- 
monics. However, one can think about an appropriate ex- 
perimental set-up, in which the Kondo couplings J"'^(i) 



oscillate with a frequency fii = ^/p, p G N. The time de- 
pendent parts in it) drive the system away from the 
isotropic two channel Kondo problem to the anisotropic 
one though their amplitudes are small compared to the 
time independent parts. Due to the oscillations in the 
input parameters, the interacting Green's functions os- 
cillate on time with a frequency of Ml/p; thus the prob- 
lem cannot be solved exactly. However, these fast oscil- 
lations allow one to average Green's functions and ob- 
servables in a period 2Tip/Vl. We have proposed a non- 
perturbative approximation, in which we have smoothed 
the fine fast oscillations around its average form. Thus, 
the Green's functions and the observables are averaged. 
Our averaging method gives the results, which are closed 
to the exact one. The higher the frequency, the bet- 
ter the approximation result. It works well for the cases 
^ ^bTk ^ Tso ^ ^at/s/t in which we study Kondo 
satellite peaks. 

The center result of our work is the satellite peak split- 
ting. This feature happens when oscillation parts are 
added to the Kondo couplings in a particular situation. 
The distance between two peaks, which are split from a 
satellite, is 2Ml/p^ the distance between two peaks, which 
are split from two consecutive satellites, is {p — 1) hfl/p, 
while the distance between two satellites is Sfi. All 
the new distances in our results depend on the number 
p, which describes the difference between the frequency 
of the Kondo couplings and the frequency of input mi- 
crowaves. We see the two close peaks clearly at small 
magnetic field in G {H) , x {H) characteristics and small 
dc voltage in G {Vdc) , X i^dc) characteristics. 

Magnetic field and applied dc and ac voltages suppress 
the Kondo correlation alternatively. When a dc source- 
drain voltage Vdc is applied, the equilibrium Kondo peak 
is split into two peaks 01 • The distance from the orig- 
inal peak to one of the new ones is dpeak = eVdc- When a 
magnetic field B is applied, the Kondopeak is also split 



into two peaks with dp 



When an 



ac source-drain voltage is applied, the Kondo peak is ex- 
panded into satellites with When 
all above external fields are applied to a SET, the equi- 
librium Kondo peak is split into satellites. The position 
of a satellite is ±eVdc ± IJ-BgiB ± (n ± 1/p) Ml when the 
Kondo couplings are considered oscillating on time with 
frequency Ml/ p. 

Although some works, in which the time dependent ex- 
change interactions like Kondo couplings J (t) were con- 
sidered, have been done by solving flow equations (3l] |. 
the exact solution for time dependent Kondo model still 
remains. In our work, we proposed the averaging method, 
which is good for high frequency regime. We predict that 
another method is needed for the adiabatic regime. Be- 
sides, one can discuss more carefully about the effect of 
microwave irradiation on the Kondo channels. Moreover, 
our work suggests qualitatively experimentalists set up 
experiments in which more than one frequency is acti- 
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vated [13. 
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